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Coherent resonant tunneling through an array of quantum dots in an inhomogeneous magnetic 
£^ ' field is investigated using an extended Hubbard model. Both the multiterminal conductance of 

an array of quantum dots and the persistent current of a quantum dot molecule embedded in 
an Aharanov-Bohm ring are calculated. The conductance and persistent current are calculated 
analytically for the case of a double quantum dot and numerically for larger arrays using a multi- 
terminal Breit-Wigner type formula, which allows for the explicit inclusion of inelastic processes. 
Cotunneling corrections to the persistent current are also investigated, and it is shown that the sign 
of the persistent current on resonance may be used to determine the spin quantum numbers of the 
ground state and low-lying excitated states of an artificial molecule. An inhomogeneous magnetic 
field is found to strongly suppress transport due to pinning of the spin-density-wave ground state 
of the system, and giant magnetoresistance is predicted to result from the ferromagnetic transition 
• induced by a uniform external magnetic field. 



PACS numbers: 73.20.Dx, 73.40. Gk, 75.30.Kz, 85.30.Vw 



I. INTRODUCTION 

Interest in the problem of coherent resonant tunneling through an. interacting mesoscopic system has been stimu- 
O | lated by a series of elegant Aharanov-Bohm (AB) ring experiments,!!! which measured the phase of the transmission 
amplitude through a quantum dot in the Coulomb blockade regime. Several theoretical works have addressed the role 
y—i \ of phase coherence in resonant tunneling through a single quantum dot: both the conductanceo and the persistent 
^ i currentu of a quantum dot embedded in an AB ring have been calculated. Many features of the experiments of 
t-H ' Ref. [j] have been explained by these model calculationsEm; however, the correlations observed between the phases of 
conductance resonances rather widely separated in energy do not appear to be explicable within these simple models, 
and it is therefore of interest to investigate coherent resonant tunneling through complex mesoscopic systems with 
, nontrivial substructure. 

In this article, we investigate the multiterminal conductance of an artificial molecule of coupled quantum dote aawell 
| as the persistent current of an artificial molecule embedded in an AB ring. Arrays of coupled quantum dotsui!3 can 
be thought of as systems of artificial atoms separated by tunable tunnel barriers. The competition between intradot 
charge quantization effects, or Coulomb blockadejLJ due to the ultrasmall capacitance of each quantum dot to its 
environment, aad-coherent interdot tunneling has been predicted to lead to a rich spectrum of many-body effects in 
these systems. E3ll3 For example, the Coulomb blockade of the individual quantum dots was predicted to be destroyed 
completely when interdot tunneling exceeds a critical value.EHl This phenomenon can be considered a finite-size anaL 
of the Mott-Hubbard metal-insulator transitionj!3 and has been observed experimentally in double quantum dotsE' 
\ A detailed theoretical investigation of the double quantum dot in the limit of a continuous energy spectrum on each 
dot, nvalid for very large dots, was able to describe the crossover from two isolated dots to a single, larger dot quite 
well.E-3 An alternative approach using an effective Hubbard-like model to describe the low-lying electronic states of the 
double quantum dot, appropriate in the limit of a discrete energy spectrum on each dot, was_also able to reproduce 
the experimentally observed crossover, as well as the nonlinear conductance of the system.li3 Here we will employ 
such an effective Hubbard model to describe coherent resonant tunneling through a one-dimensional (ID) array of 
quantum dots coupled to multiple leads. We will briefly consider the corresponding situation for the coherent resonant 
tunneling through a two-dimensional (2D) quantum dot array as well. We focus on the strongly-correlated regime, 
where interdot tunneling is too weak to destroy the energy gap stemming from Coulomb blockade effects, and where 
the energy spectrum on each dot is discrete. 

An important consequence p£ coherent interdot tunneling in the strongly-correlated regime is the formation of 
interdot spin-spin correlations!!! analogous to those in a chemical bond at an energy scale J ~ t 2 /U, where U 
is the charging energy of a quantum dot and t — (Ti 2 /2m*) J d 3 x 5 , ^ l (x)V 2v I / n (x) is the interdot hopping matrix 
element, ^ m ,n being electronic orbitals on nearest-neighbor dots. In a system with magnetic disorder, such a spin 
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configuration is pinned, and the resulting blockage of spin backflowlIS leads to strong charge localization. However, an 
applied magnetic field will break such an antiferromagnetic bond when the Zeeman splitting gfisB > </, leading to 
an enormous enhancement of the charge mobility. Such spin-dependent many-body effects on the magnetotransport 
should be experimentally observable provided r+T'*), ksT < J, where r + T^ is the total broadening of the resonant 
levels of the system due to finite lifetime effects and inelastic scattering; they can be readily distinguished from 
orbital effects in arrays of quasi-two-dimensional quantum dots by applying the magnetic field in the plane of the 
dots. Observation of the predicted giant magnetoresistance effect in the low-temperature transport through coupled 
quantum dots would, we believe, represent a clear signature of the formation of an artificial molecular bond. 

The paper is organized as follows: In Sec. II, an extended Hubbard model describing the low-lying electronic states 
of an array of coupled quantum dots is introduced, and the magnetic phase diagram of the system is discussed. In 
Sec. Ill, some general expressions for the conductance and persistent current of an interacting mesoscopic system 
coupled to multiple electron reservoirs are derived. In Sec. IV, the conductance through a double quantum dot in an 
inhomogeneous magnetic field is calculated. The magnetoresistance of the system is shown to be proportional to U 3 /t A . 
The persistent current through a double quantum dot embedded in an AB ring is also investigated. For a system 
with an odd number of electrons, it is shown that resonant tunneling through molecular states with odd and even 
S z leads to contributions to the persistent current of opposite signs. In Sec. V, the conductance of one-dimensional 
arrays of quantum dots is investigated. The inelastic scattering rate in the system is shown to be suppressed in large 
arrays due to the orthogonality catastrophe. A pronounced suppression of certain resonant conductance peaks in an 
applied magnetic field is predicted to result from a field-induced ferromagnetic transition. A many-body enhancement 
of localization is predicted to give rise to a giant magnetoresistance effect in systems with spin-dependent disorder. 
Such spin-dependent magnetoresistance effects are found to be much weaker in the ballistic transport regime. In Sec. 
VI the conductance in the 2D quantum dot array is investigated, and the qualitative physics is found to be similar to 
that in the ID array. Some conclusions are given in Sec. VII. 



II. MODEL 

The system under consideration (Fig. [I]) consists of a linear array of quantum dots electrostatically defined! - ! in a 
2D electron gas, coupled weakly to several macroscopic electron reservoirs, with a magnetic field in the plane of the 
dots. Each quantum dot is modeled by a single spin-1/2 orbital, representing the electronic state nearest the Fermi 
energy Ef, and is coupled via tunneling to its neighbor and to one or more electron reservoirs. Transport .occurs 
between the left (L) and right (R) reservoirs; reserjioirs 1 to Nd are considered to be ideal voltage probesjijpajjui 
serve to introduce inelastic processes in the system.E3 Electron-electron interactions in the array are describecO'EH 
by a matrix of capacitances Cjj : We assume a capacitance C g between each quantum dot and the system of metallic 
gates held at voltage V g , an interdot capacitance C, and a capacitance C r between a quantum dot and each of its 
associated electron reservoirs. The diagonal elements of Cy are the sum of all capacitances associated with a quantum 
dot, C\\ = CN d N d = C g + C + 2C r , Ca = C g + 2C + C r (1 < i < Nd) and the off-diagonal elements are CV,- = — C 
for nearest neighbor dois.^'. These capacitance coefficients may differ from their geometrical values due to quantum 
mechanical corrections ,oEa but enter only as parameters in our model. The Hamiltonian of the quantum dot array is 

H dots = ^d) a d ja + (ti«4 +la d ja + H.c.) + 1 Y^(Qi + Q g )Cr.\Q j +Qg )-9lj2 C%\ (1) 

i,<* o," *,j »,j 

where d^ a creates an electron of spin a in the jth dot, Qj = —e dj a dj a is the charge operator for dot j, Q g = C g V g 
is a polarization charge induced by the gate, and 

e jo • mi, ■> <r,^V ; ,'i;. (2) 

i 

where Bj is the Zeeman splitting on dot j and ej is the orbital energy of the quantum-confined orbital under consid- 
eration on dot j. The last term in Eq. (||) represents a shift in the orbital energy due to the capacitive coupling to 
the reservoirs. This term has an important effect on the nonlinear transport,E3 but does not affect the linear response 
and equilibrium properties which aije-ithe subject of the present article. We therefore set C r = in the following. 

-ffdots reduces to a Hubbard modefca with on-site repulsion U — e 2 jC g in the limit C — > 0. In general, -f/dots describes 
an extended Hubbard model with screened long-range interactions. The elements of the inverse capacitance matrix 
decrease exponentially with a screening length that increases with C/C g . For C <C C g , C~^ ~ {C/C g y % ~^ /C g , while 

for C/Cg — > oo, C^j 1 — > 1/NdCg, and intradot charging effects are fully screened. Electronic correlation effects are 
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thus decreasing functions of C/C g . By varying the interdot electrostatic coupling, one.-can thus study the transition 
from a strongly correlated artificial molecule exhibiting collective Coulomb blockadeO for C -C C g to a ballistic 
nanostructure where correlation effects are negligible for C>C 9 . 

Of particular interest to us here is the magnetic phase diagram of the quantum dot array. In the strongly-correlated 
regime, where intradot charging is not strongly screened, the TV-electron ground state of i?dots will form a spin-density- 
wave (SDW) due to interdot superexchange.O In an external magnetic field, the electron spins will tend to align with 
the field to minimize the Zeeman energy. There is thus a transition from a spin-density-wave ground state at B = 
to a -ferromagnetic state at some critical magnetic field B c . For an infinite ID array with C = and ( « (7, one 
findsEfl 

At 2 

g^BB c ~ — -(27m - sin 27m), (3) 
iru 

where n = 1 — |1 — N/L\ is the filling fraction of electrons (holes) for TV < L (N > L) in the orbital under 
consideration. Recall that we are here considering only the single spin-1/2 orbitaLnearest Ep in each quantum dot; 
the magnetic field required to spin-polarize an entire quantum dot is much larger £3 In a system with C > 0, intradot 
Coulomb interactions are screened, and B c is therefore expected to increase. Fig. || shows the spin susceptibility Xs 
for C/Cn — 1/2 in linear arrays with 8 electrons on 12 dots and 10 electrons on 10 dots. The n-dependence of B c 
in Fig. @ is qualitatively similar to that in a system with intradot interactions only, but the values of B c are roughly 
twice those of a system with C = 0. Note the rapid growth of Xs as B — > B c . In an infinite,-array, Xs is expected to 
diverge as B — > B c because the system undergoes a second order quantum phase transition.^ The spin-polarization 
transition (SET) predicted to occur in an array of coupled quantum dots is in contrast to that observed in a single 
quantum dot,E3 where the critical point occurs for minimum total spin. 

In the following, we shall investigate the effects of the SDW correlations and the SPT on the low-temperature 
magnetotransport through coupled quantum dots. 



III. QUANTUM TRANSPORT FORMALISM 

Before investigating the particular Hamiltonian of interest, Eq. (Q), it will first be useful to derive some general 
formulas for the conductance and persistent current arising due to resonant tunneling through an interacting system. 
An arbitrary interacting mesoscopic conductor coupled to M macroscopic electron reservoirs is described by the 
Hamiltonian 

M M 

h = H int ({di d n }) + ekc i ck + E E E ( v ^Ad n + h.c) , (4) 



a=lk£a a=l k£a n 



where creates a complete set of single-particle states in the mesoscopic system, c| e creates an electron in state 
k of reservoir a, and H[ nt is a polynomial in {d\, d n } which commutes with the electron number TV = ^ n d n . Here 
the spin index a has been absorbed into the subscripts n and k. We denote the ground state of -Hint for each N by 
| Oat) and the ground state energy by E^. We assume E^ to be nondegenerate, as is generically the case in a nonzero 
magnetic field. 

If the tunneling barriers to the reservoirs are sufficiently large, and if the temperature and bias are small compared 
to the energy of an excitation, then the main effect of particle exchange with the reservoirs will be to cause transitions 
|0jv-i) — > |0jv) between the nondegenerate ground states of the system. In the vicinity of such a resonance, the 
nonequilibrium Green's functions describing [propagation within the system in the presence of coupling to the leads 
can be shown to have the Breit-Wigner formEj 

i \ (0/Y-lMn|0/Y)(0/Y|d^|0/v_i) ..... . . 

G r nm (0 = " ' o ' P0 — T r /o + additional poles, 5 

e - E" N + E^_ 1 + iT N /2 



(0 N -i\d n \0 N )(0 N \dl\Q N -i) E Q Wa(e) 



G wn( e ) = 7 F o , F o m « i /r %\2 + additional poles, (6) 

[e-E N + E N _ 1 )- + (1 N /2) 

where f a (e) = {exp[(e — /i Q )/fc/jT] + 1} _1 is the Fermi function for reservoir a, Tjv = Y] a — 1 and 

- 27rJ2T / ( ^- 1 \ V ^\° N )i° N \ V krndU0N-i)S(e k -E% + (7) 
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Here, G^£[(e) are Fourier transforms of the Keldysh Green's function Gf im {t) — i{d\ ri {Q)d n {t)) and the retarded Green's 
function G r nm {t) = — iG(t){{dn(t), d> m (Q)}), respectively. With the aid of these Green's functions, the conductance and 
persistent current resulting from resonant tunneling through the system can be calculated. 



A. Multiterminal Conductance Formula 

The expectation value of the currentJiowing out of the interacting region into reservoir a can be expressed using 
the formalism of Meir and Wingreen asE3 

I * = -lj deImTr{r a (e)[G<(e) + 2f a (e)G r (e)]}, (8) 

where r" m (e) = 27r^ fcgQ Vk n V£ m 5(e — ek) is a matrix characterizing the tunnel barrier connecting reservoir a to the 
system. laserting G^£[ from Eqs. ((^) and (JsJ) into Eq. (Js[) , one finds the multiprobe current formula for resonant 
tunnelinga 

" h H ] n (e~E% + + (T N /2f 

The low-temperature transport through such a correlated many-body system weakly coupled to multiple leads thus ex- 
hibits resonances of the Brcit-Wigner type,E3 where the positions and intrinsic widths of the resonances are determined 
by the many-body states of the system. Eq. (|^), which expresses the current in terms of transmission probabilities, is 
a generalization of the multi-terminal conductance formula for a noninteracting system derived by BiittikercH to the 
case of resonant tunneling through an interacting system. 

In deriving Eq. (||), we have neglected the additional poles in G<£(e), which is justified provided fcgT, <C 
AE N and < AE N , where AE N = min(£^ - - E%_ X ,E 6 N+1 - E°- /i Q , fi a - E%_ x + E° N _ 2 ), Ejf 

being the energy of the lowest lying excited state of the iY-eiectron system. Eq. (g) is thus appropriate to describe. 



resonant tunneling through semiconductor nanostructurescSEa or ultrasmall metallic/superconducting systems^Hil 
under conditions of low temperature and bias, where transport is dominated by a single ground-state to ground-state 
transition |0jv-i) <-> \0n)- Eq. (^) is not applicable to systems with a (spin)degenerate ground state [AE(N) = 0], 
where the low-temperature physics is that of the Kondo effect, as discussed in Refs. |2l^j32| -|34|. 

We next specialize to the configuration shown in Fig. |l|. Transport occurs between the left (L) and right (R) 
reservoirs. The auxiliary reservoirs 1, . . . , Nd are assumed to be connected to ideal voltmeters .Efl An ideal voltmeter 
should have an infinite impedance, so we demand that the expectation value of the current flowing into reservoirs 
1, . . . , Nd be zero, which fixes n\, . . . , fj,N d via Eq. (|9|). Eliminating /i(e), . . . , fw d (e) from Eq. (|^), and taking the 
linear response limit, one finds the effective two-terminal conductance between the left and right contacts 



', r i't r » J (e - E% + E°,_i) 2 + (I>/2) 2 



The total width of the iVth resonance may be written Tjv = + Tfj + ' , where = Yla=i ^n- The quantity 
Tj9 /% may-, be interpreted as the total inelastic scattering rate due to phase-breaking processes in the auxiliary 



Such processes arise when an electron in the ith dot escapes into reservoir i and is replaced by another 
electron from the reservoir, whose phase is uncorrelated with that of the previous electronic 

For simplicity, let us assume that the tunnel barriers coupling the system to the external reservoirs are described by 
the energy-independent parameters 27r^] fceQ , |Vfc„| 2 £(e fc — E) = T^S na , a = 1, . . ., N d ; T5 nl , a = L; TS nNd , a = R. 

Then the partial widths of the Nth resonance are simply T% = TJ2 a I(0a'MLI a'-i)| 2 > t n = r T,* I (°N Mjv d J JV-i) P> 
and 

N d 

r% = t^Y,WM^n-i)\ 2 . (ii) 

i=l a 

Since {dl a } creates a complete basis of single-particle states in the array, one would have = for a noninteracting 
system. However, in an interacting system the wavefunction ovsriap |(0jv|<if o .|0jv-i)| is suppressed by correlation 
effects, leading to an orthogonality catastrophe in the large- JV limit.Ea The many-body suppression of the wavefunction 
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overlaps leads to a reduction of both the elastic broadening Ty = + and of the inelastic broadening Ty/ of 
the conductance resonances, so that the system becomes more and more weakly coupled to the environment as N 
increases. The suppression of the wavefunction overlap with increasing N is shown for an array of ten quantum dots 
with pure Hubbard interactions in Fig. ^. The effect of-such an orthogonality catastrophe in the seqoiential tunneling 
regime has previously been discussed by Kinaret et al£3 and by Matveev, Glazman, and BarangerO 



B. Persistent current 



Let us next consider the persistent current through the quantum dot array when the left and right reservoirs in Fig. 
1(a) are connected together to form a ID ring of length L enclosing an AB flux $ = (Tic/e)(f>. Since the persistent 
current is an equilibrium property, the auxiliary reservoirs must be in mutual equilibrium at electrochemical potential 
fx, as indicated in Fig. 1(b). The persistent current in such an open system is given by /(</>) = —{e/Ti)dVL/ d<p, where 
f2 is the grand canonical potential. The grand canonical potential-may be determined from the electronic scattering 
matrix of the system, as discussed by Dashen, Ma, and BernsteirO and by Akkermans et al.^B and is 

n W = _ /^/(£)Im In det £(£,<£), (12) 

where S(E, </>) is the scattering matrix of the multiply-connected structure shown in Fig. 1(b) and f(E) is the Fermi 
function with electrochemical potential /z. In order to obtain S{E,<p), it is useful first to introduce the scattering 
matrix S(E) for the structure shown in Fig. 1(a), which is related simply to the retarded Green's function, 

S W (E) = -6 kk , + 2m %^G; m (^), (13) 

where Vkea is the group, velocity of state k in reservoir a. Using Eq. (§), one sees that S(E) has the Breit-Wigner form 
discussed by ButtikerEJ S(E) may be divided into submatrices rkk' , describing reflection of mode k' in the reservoir 
back into mode k in the reservoir, e k and 5k, describing transmission of mode k from the reservoir into the the ring to 
the left and right, respectively, 7, describing transmission of a circulating state in the ring through the quantum dot 
array, and a and (3, describing reflection of the circulating states in the ring at the left and right ends of the quantum 
dot array, respectively. In terms of the submatrices of S(E), the flux-dependent scattering matrix for the combined 
structure may be written 

c, tTP ,s . a5 k S k ' + f3e k s k > - (7 - cos</>e~ 4pi )(e fc 4' + S k e k ') 

S kk > {E, 4>) = r k k' H 2 a — o 1 -ivL 1 — =2»I ' 14 

7^ — ap — Z7 cos tpe % v L + e Zi p l 



where p — \ / 2mE/U and it is assumed that ek = e k ' = E. Substituting Eq. (|14| ) into Eq. (|i2|), and taking the derivative 
with respect to <j), one obtains the persistent current due to coherent resonant tunneling through the quantum dot 
array. The general expression for I((f>) is somewhat cumbersome, but for the special case = T^j = T^/2, one 
finds the simple result 



r(e) . dEf(E)T% \ ^J^P L 

eTy f (cos0 - cospL) 2 



2h 



_ (r^ ) /2)sinpL 
N N — l cospL — cos 



E - E° + E%_, - 



T^/2) 2 



(15) 



The persistent current is thus also determined by the conductance matrix elements T^. From Eq. ([lfj) one sees that 
as the energy of the resonance is tuned through the Fermi level (by varying the gate voltage V g ), the persistent current 

exhibits a maximum of height ~ (eT$/h) sine/) sin kpL, which decays with a width Tn when the resonance is moved 
above /1. / also decays when the resonance is moved below the Fermi level, but in an oscillatory pattern determined 
by the level spacing in the ring. 

In the limit (closed system), the integrand in Eq. ( |l5| ) consists of a series of delta functions whose weights 

give the current contributed by each of the discrete states in the ring which couple to the dot array, subject to periodic 
boundary conditions. If the ring has dimensions comparable to those of the dot array, then the level spacing in the 
ring may exceed the width of the resonance (Kvp/L > Tjv), and the above approach will break down. The persistent 
current will then be dominated by charge fluctuations coupling the highest occupied state in the ring, of energy ep, 
with the lowest unoccupied many-body state in the dotarrayB The coupling matrix element is 
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\t\ 2 =t% + t 2 L +2t R t L cos <p, (16) 

where the matrix elements 

t L =J2(°N\V kF id{ lT \0N-i) (17) 

(J 

tR = E^I^JV^JOjV-!) (18) 
er 

may be chosen real. A straightforward calculation! then yields the persistent current through the quantum dot array 

T( j.\_ 2(e/h)t R t L sin0 

W [(tF-E% +n-i) 2 +i\t\ 2 ] 1/2 ' [ ' 

The product tpti, may be greater or less than zero depending on the relative signs of Vk F i and Vk F N d , and on the 
relative signs of the wavefunction overlaps, so the persistent current on resonance may be either diamagnetic or 
paramagnetic. For a purely ID system of spinless electrons, or with an even number of spin-1/2 electrons, the sign of 
I(4>) is determined by the total, number of electrons in the system, as discussed in Rcf. This is a manifestation of 
the so-called Leggett theoremEj However, in general the sign of I(cf>) depends on the particular state which dominates 
the resonance, and may be used to classify the quantum numbers of that state, as discussed below. 



IV. DOUBLE QUANTUM DOT 



Let us first consider the case of double quantum dot, for which the conductance matrix elements (0) can be obtained 
analytically. For Nj = 2, i/dots reduces to a two-site Hubbard model with on-site repulsion U — e 2 Cs/(C|. — C 2 ) 
and nearest neighbor repulsion Uu = e 2 C/(C 2 : — C 2 ), where Cs = C g + C + 2C r - Spin disorder is introduced via 
a Zeeman splitting on dot 1, B\ = At A, B 2 = 0. Experimentally, such an inhomogeneous field could be produced, 
e.g., by the presence of a small ferromagnetic particle. The total width of the one-particle resonance is found to be 
ri = T + rW, and the prefactor in Eq. @ is rfTf /(Tf + Tf ) = (T/4)/(l + A 2 ) = T . The maximum conductance 



■ rf>- 

at the one-particle resonance is thus 



e 2 /h 



T = 



(i + A 2 )(i + r«/r)' 
e 2 r /Ahk B T, r + r (i) < k B T < Wi + a 2 . 



(20) 



Inelastic scattering suppresses the resonant conductance at T = 0, but has no effect when the resonance is thermally 
broadened. For U — XJ\i S> t, the two-particle ground state of the double quantum dot has an antiferromagnetic spin 
configuration characterized by the superexchange parameter 



J = 2thy - A + ^7 2 + A 2 j , (21) 

where 7 = t/(U — U\2)- Note that 2t^j < J < 4^7. The two-particle resonance is separated from the one-particle 
resonance by eAQ g /(Cs — C) = U12 + 2t(l + A 2 ) 1 / 2 — J, and the conductance is determined by the matrix elements 



(02|dJ T |0i) 
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At VA^ + l) 



1/2 



1 T A/ 7+ ^7 2 + A 2 



1 + A ± Va^TT) 



1/2 



(22) 



where A 2 — 1 + A 2 / (7+ -\/7 2 + A 2 ) 2 and the upper (lower) sign holds for j = 1 (2). For B\ > J, the antiferromagnetic 
spin configuration is pinned, leading to a strong suppression of the amplitude to inject electron 2 into dot 1, and a 
concommitant suppression of the second conductance peak. Inserting Eq. ( |22| ) into Eqs. (Q) and (|l0|), one finds the 
T = resonant conductance 



g; = i6( e 2 // l )( 7 /A) 2 /(i + r«/r), i « a « 1 
= 4(e 2 //i) 7 2 /(i + r( l Vr), A»l. 



(23) 
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A second doublet of conductance peaks for N = 3, 4 is separated from this doublet by AQ g ~ e (center to center), 
and one finds G3 = G 2 , G\ = G\ due to electron-hole symmetry. The resonant conductance for N = 2 is suppressed 
by a factor of j 2 compared to that for N — 1 due to collective spin pinning (one readily verifies that the resonant 
conductance is suppressed by the same many-body factor in the regime of thermally broadened resonances). This 
dramatic many-body suppression of the conductance is illustrated in Fig. ^ for several values of 7. The effect of spin 
disorder is to be contrasted with that of a charge detuning A = (ei — £2)/2t, investigated by Klimeck et alS3 and by 
van der Vaart et al.p for which both G\ and G 2 are given by Eq. ( |20| ) at T = (G 2 is then reduced by a factor of 
2 in the thermally broadened regime) . The very different effects of spin and charge disorder stem from the fact that 
the repulsive interactions in Eq. (Q) enhance spin-density fluctuations, but suppress charge-density fluctuations. 

Let us now consider the effect of an additional homogeneous magnetic field applied parallel to the inhomogeneous 
field, Bi = 4tA + B, B 2 = B. For B > J, it is energetically favorable to break the antiferromagnetic bond between 
the dots and form a spin-polarized state, thus preventing collective spin pinning effects. G 2 is then given by Eq. (pc|). 
The resulting magnetoresistance on resonance for T = and A ^> f is thus 

AR* _ h gn B (i) h g^ B {U - t/12) 3 (0A , 

AB ~ e 2 4J 7 2l + e 2 8i 4 ' [ ' 

In the thermally broadened resonance regime, the factor (I + /T) is replaced by 2k B T/nT. Since the Coulomb 
energy U — U\ 2 is typically large compared to the interdot tunneling matrix element t, the predicted magnetoresistance 
is extremely large. This giant magnetoresistance effect is a direct indication of the field-induced breaking of the 
artificial molecular bond between the dots.Ell 

The conditions necessary to observe the predicted magnetoresistance effect may be determined by including the 
effect of transport through the triplet excited state via the method of Refs. [ll] and [l3|. One finds the resonant 
conductance at B = for k B T > T + 

= e 2 exp(/3J) / 2T t \ 

2 2hk B T 2 exp(/3 J) - I V exp(/3J) + I J ' K ' 

where T 8 ~ 4 7 2 T is the sequential tunneling rate through the pinned antiferromagnetic ground state and T t = Tq is 
the sequential tunneling rate through the triplet excited state. The magnetoresistance is thus reduced by a factor of 2 
at a temperature k B Ti/ 2 = J/ ln(r t /T s ). Increased coupling to the leads and/or inelastic scattering can be shown to 
lead to a similar admixture of transport through excited states when r + F 1 ' ~ J. We therefore expect the predicted 
giant magnetoresistance effect to be observable for k B T, F + rW < J. In currently available GaAs quantum dot 
systems, charging energies are typically of order lmeV, and one expects tunneling matrix elements t ~ O.lmeV for 
moderate to strong interdot tunneling, so values of J in the range .01 — O.lmeV should be attainable. 

Let us next consider the persistent current through the double quantum dot. From Eqs. (|T^) and (|l9|), one sees 
that the persistent current is also suppressed at the N = 2 resonance due to the many-body factor (|22|). However, 
in the later case of a nanoscopic ring with level spacing Kvf/L 2> the suppression of the persistent current is 
only linear in (02|cZL|0i). An interesting question is the effect of cotunneling through excited states of the double 
dot when the dot ring-coupling \t\ exceeds the many-body level spacing in the double dot J. In order to address this 
question, we have studied the closed system of a double quantum dot embedded in an AB ring numerically using the 
Lanczos techniquecll (see Fig. |[). The Hilbert space was truncated by discretizing the ring (8 sites were used). In 
order to distinguish the contributions from tunneling through the S z = ground state and the S z = 1 excited state 
of the double dot, the total number of electrons in the system was chosen to be odd (in this case 5). Thus if the tatal 
number of up-spin electrons iVj is even, the total number of down-spin electrons N± must be odd, and vice versa.Ej In 
the weak-coupling limit \t\ <C J, where a single level of the double dot contributes to the resonant current, the spin 
of the tunneling electron is well defined, and is 

ha/2 = (0 N \S Z \0 N ) - (Ow-il^lOjv-i). (26) 

Fig. H shows the persistent current at (f> = ir/2 as a function of the gate voltage Q g in the vicinity of the first Coulomb 
blockade doublet centered near Q g = e/2. The doublet splitting is here enhanced due to the finite level spacing in 
the ring. For A > 0, the ground state of the coupled dot-ring system generally has 2Vj_ = 2Vf + 1 (in this case Ni = 3 
and -ZVj = 2). The first electron to enter the double dot as Q g is increased from zero enters the lowest single-electron 
eigenstate o£-jbhe double dot, and thus has a = — 1. Since N± is odd, the resonant current is diamagnetic due to the 
parity effectfl The second electron to enter thexlouble dot goes into the state IO2) and thus has a = +1. Since N-\ is 
even, the resonant current is thus paramagnetica (see solid curve in Fig. |^) . The height of the second peak is reduced 
compared to that of the first, but by a smaller factor than for the conductance [c.f., Eqs. ( p^| ) and (|Tjj)]. However, 
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there is also a contribution to the persistent current due to cotunneling through the first excited state of the double 
dot, which is higher in energy by J than | O2} ■ This state has a = — 1, and thus leads to a diamagnetic contribution to 
/(</>). This state couples more strongly to the leads, but is suppressed by a large energy denominator when \t\ <C J- 
As \t\ is increased (dotted and dashed curves in Fig. ^), the cotunneling contribution becomes increasingly important, 
and there is a crossover from a small paramagnetic peak to a larger diamagnetic peak for \t\ > J. Fig. a clearly shows 
that the sign of the persistent current induced by tunneling through a ID structure may be used to characterize the 
spin quantum numbers of the ground state and low-lying excited states of such a system. 

V. ID ARRAY OF QUANTUM DOTS 

Let us next consider tunneling through largeij-arrays of quantum dots. For AT d > 2, the TV-body ground states of 
Eq. (|l|) were obtained by the Lanczos technique,^ and the conductance was calculated using Eq. (|l^). At T = and 
in the absence of inelastic scattering, the conductance peaks all have height e 2 /h in the absence of disorder, since 
in that case F^ = F^ = Fjv/2. Inelastic scattering leads to additional broadening of the conductance peaks, and 
suppression of the resonant conductance below e 2 /h. Disorder also leads to a suppression of the T = resonant 
conductance below e 2 /h due to the breaking of left-right symmetry ^ T 1 ^. In the following, we concentrate on the 
thermally broadened resonance regime, where the peak heights of the conductance resonances depend most strongly 
on the conductance matrix elements and T^. For UbT Tn, Eq. ( |Toj ) simplifies tcErES 

G = T E \-fb -z° N + JSjU)]. (27) 

'* N 1 N ~r 1 N 

Fig. |FJ shows the conductance through a linear array of 10 quantum dots with C = as a function of the chemical 
potential \x in the leads, whose value relative to the energy of the array is controlled by the gate voltages. The two 
Coulomb blockade peaks in Fig. || are split into multiplets of 10 by interdot tunneling, as discussed in Refs. 
refer to these multiplets as Hubbard minibands. The energy gap between multiplets is caused by collective Coulomb 
blockade,E3 and is analogous to the energy gap in a Mott insulator a The heights of the resonant conductance peaks 
in Fig. |6|(a) can be understood as follows: Since the barriers to the leads are assumed to be large, the single-particle 
wavefunctions of the array are like those of a particle in a one-dimensional box. The lowest eigenstate has a maximum 
in the center of the array and a long wavelength, hence a small amplitude on the end dots, leading to a suppression of 
the 1st conductance peak. Higher energy single-particle states have shorter wavelengths, and hence larger amplitudes 
on the end dots, leading to conductance peaks of increasing height. The suppression of the conductance peaks at the 
top of the 1st miniband can be understood by an analogous argument in terms of many-body eigenstates; the 10th 
electron which enters the array can be thought of as filling a single hole in a Mott insulator, etc. 

In Fig. ||(b), the spin-degeneracy of the quantum dot orbitals is lifted by the Zeeman splitting. There is a critical 
field B c above which the system is spin-polarized [c.f., Eq. (||)] Because B c is a function of n, one can pass through 
this spin-polarization transition (SPT) by varying n at fixed B. In Fig. ^|(b), this transition occurs between the 4th 
and 5th electrons added to the array, consistent with the prediction of Eq. (0) . The effect of this transition on the 
conductance spectrum is dramatic: The first 4 electrons which enter the array have spin aligned with B (up), but 
the 5th electron enters with the opposite spin, and goes predominantly into the lowest single-particle eigenstate for 
down-spin electrons, which couples only weakly to the leads, leading to a suppression of the 5th resonant conductance 
peak by over an order of magnitude. It should be emphasized that the heights of the conductance peaks change 
dis ■continuously as a function of B each time there is a spin-flip. 

Splitting of the Coulomb blockade peaks due to interdot coupling and suppression of the conductance peaks .at the 
miniband edges have recently been observed experimentally by Waugh et alu However, it has been pointed outl3 that 
both effects can also be accounted for by a model 2 -!! of capacitively coupled dots with completely incoherent interdot 
tunneling. It is therefore of interest to consider the effects of interdot capacitive coupling in the regime of coherent 
interdot transport. A nonzero interdot capacitance C introduces long-range electron-electron interactions in Eq. (Q) 
and decreases the intradot charging energy U. Fig. || shows the spin susceptibility \s for CJC g = 1/2 in linear arrays 
with 8 electrons on 12 dots and 10 electrons on 10 dots. The n-dependence of B c in Fig. [2] is qualitatively similar to 
that in a system with intradot interactions only, but the values of B c are roughly twice those of a system with C — 0. 
Note the rapid growth of \s as B — > B c . In an infinite-Array, Xs is expected to diverge as B — > B c because the system 
undergoes a second order quantum phase transition.Ej ThepSPT predicted to occur in an array of coupled quantum 
dots is in contrast to that observed in a single quantum dotp3 where the critical point occurs for minimum total spin. 

Disorder introduces a length scale which cuts off the critical behavior as B — > B c . However, as shown in Fig. 
f^, where disorder St ~ t has been included in the hopping matrix elements, the SPT has a clear signature in the 
magnetotransport even in a strongly disordered system. In Fig. [?], the peak splitting due to capacitive coupling is 
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roughly ten times that due to interdot tunneling, so that the peak positions are within ~ 10% of those predicted by a 
classical charging modeLEa However, the dramatic dependence of peak heights on magnetic field — the 4th conductance 
peak in Fig. 0(b) is suppressed by a factor of 32 compared to its B = value due to the density-dependent SPT 
described above — can not be accounted for in a model which neglects coherent interdot tunneling. This effect should 
be observable provided g/j,BB c > max(fcsT, Ti/r,), where is the inelastic scattering time. We believe that this 
striking magnetotransport effect is the clearest possible signature of a coherent molecular wavefunction in an array of 
quantum dots. 

Fig. |^ shows the conductance spectrum for an array of 6 quantum dots with the same parameters as in Fig. 0, but 
with spin-dependent disorder in the hopping matrix elements, as could be introduced by magnetic impurities. Several 
conductance peaks at B = (solid curve) are strongly suppressed due to a many-body enhancement of localization. 
This effect arises becausq-*epulsive on-site interactions enhance spin-density wave correlations, which are pinned by the 
spin-dependent disorders At B = 1.3T (dotted curve) the system is above B c and is spin-polarized, circumventing 
this effect. The second conductance peak is enhanced by a factor of 1600 at 1.3T compared to its size at B = (not 
visible on this scale). This giant magnetoconductance effect is a many-body effect intrinsic to the regime of coherent 
interdot transport. 

Another interesting phenomenon stemming from the competition between coherent interdot tunneling and charging 
effects is the Mott-Hubbard metal-insulator transition (MH-MIT), which occurs when collective Coulomb blockade 
(CCB)li3 is destroyed due to strong interdot coupling. For GaAs quantum dots larger than about lOOnm in diameter, 
we find that this transition is caused by the divergence of the effective interdot capacitance, similar to the breakdown 
of Coulomb blockade in a single quantum dot £3 Within the framework of the scaling theory of the MH-MIT ,£2 one 
expects a crossover from CCB to ballistic transport in a finite array of quantum dots when the correlation length £ in 
the CCB phase significantly exceeds the linear dimension L of the array. Fig. |^ shows the conductance spectrum for 
5 quantum dots with 5 spin-1/2 orbitals per dot. The divergence of the effective interdot capacitance as the interdot 
barriers become transparent is simulated by setting C^/C g = 2" -1 , n — 0, ... ,4. In Fig. || minibands arising from 
each orbital are split symmetrically into multiplets of 5 peaks by CCB, with the center to center spacing between 
multiplets equal to e 2 /C g , while the energy gap between minibands corresponds to the band gap ~ A enhanced by 
charging effects. The CCB energy gap is evident in the first 3 minibands, but is not resolvable for the higher orbitals 
{C/Cg > 4), although there is still a slight suppression of the conductance peaks near the center of the 4th miniband. 
Comparison of the compressibility of the system to a universal scaling function for the MH-MIT calculated by the 
method of Ref. ^3| indicates £/L ~ 10 3 for C/C g — 8, so that the transport in the 5th miniband is effectively ballistic. 
The peak spacing within a miniband saturates at e 2 /LC g (plus quantum corrections ~ t/L) in the ballistic phase 
because the array behaves like one large capacitor, as observed experimentally in Ref. 0. Fig. ||(b) shows the effects 
of a magnetic field on the conductance spectrum: a sequence of SPTs is evident in the different minibands, with B c 
an increasing function of C/C g , leading to quenching of magnetoconductance effects in the ballistic regime. 

A finite-size scaling analysis of the compressibility indicates that the MH-MIT probably occurs at C ' jCg^, oo in 
an infinite array of quantum dots, when the interdot barriers become transparent to one transmission mode.Ea 



VI. 2D ARRAY OF QUANTUM DOTS 

Finally, we briefly consider coherent tunneling through a 2D quantum dot array (Fig. [To|) to investigate whether 
the many-body giant magnetoconductance effect discussed in Sec. V arises in two dimensions as well. We use Eq. fl27j ) 
to calculate the linear tunneling conductance through a 2D 3x3 quantum dot array consisting of Nd = 9 quantum 
dots. The corner dots in the array are weakly coupled to electron reservoirs as shown in Fig. [l^. The Hamiltonian 
of the array is the same as that given by Eq. (|l|), where the second term is now modified to incorporate the nearest 
neighbor tunneling in the 2D array: the tunneling amplitudes connecting two nearest neighbor dots i and j in the 

array are multiplied by the Peierls phase factors e %K h ^ dl,J with A as the magnetic vector potential.0 We consider 
a uniform flux (f> — Bo? piercing each unit cell of the array in Fig. |l^. The magnetic field B enters through the 
tunneling amplitudes iy and through the intradot single-particle field dependence entering £j CT in Eq. (|). The flux 
sensitive phase factors in Eq. (Q) lead to a flux periodic modulation of the linear conductance with periodicity given 
by one fundamental flux unit hc/e. This flux dependence of .-the linear conductance and the associated ground state 
persistent current oscillations have been discussed elsewhereO Here we follow our discussion in the previous section, 
concentrating on a fixed applied field, and focus on the magnetic field induced spin effect (i.e. the single-dot Zeeman 
physics) on the linear conductance peak heights in the 3x3 array. 

The partial width of the Nth resonance is plotted as a function of N in Fig. [ll] (a) and (b) for the applied 
magnetic field B = and B = 1.3T in the array, respectively. The corresponding linear conductance in the 3x3 array 
is shown in Fig. 121 The peak splittings in the linear conductance in the 3x3 array are not distributed uniformly, 







but the shape of the envelope function for conductance peak heights is similar to that in ID chains (e.g. compare 
Figs, [j] and [12]). This envelope function is peaked at quarter filling in the Hubbard model we use here. In the strongly 
correlated Hubbard array, the overlap matrix element is approximately given by ~ [1 — Pn-i(1)]Pn(1) with 
P/v(l) = J^OjvldiffdiirlOjv) being the probability to find the corner dot occupied. For U/t — 10, the maximum and 
minimum P/v(l) are approximately 1 and 0, leading to a peak in the linear conductance at quarter filling (and also at 
three quarter filling due to particle- hole symmetry) in Hubbard arrays. The addition spectra shown in the conductance 
vs. chemical potential plots of Figs. [ll] and [l^ correspond to a sequence of ground states that are characterized by 
the number of electrons N, the total spin S, and the component S z of the total spin along the quantization axis. 
The ground state (N,S,S Z ) sequence for B = in the array is (1,1/2,1/2) -> (2,0,0) -> (3,3/2,1/2) -> (4,1,1) -> 
(5,1/2,1/2) -> (6,0,0) -> (7,1/2,1/2) -> (8,0,0) -> (9,1/2,1/2). In the B = 1.3T array the similar sequence 
of the ground states is given by (1,1/2,1/2) -> (2,1,1) -> (3,3/2,3/2) -> (4,2,2) -> (5,5/2,5/2) -> (6,2,2) -> 
(7, 1/2, 1/2) — > (8, 1, 1) — * (9, 1/2, 1/2). It can be seen from the latter sequence that the 6th linear conductance peak 
at B = 1.3 is suppressed (by a factor of 23) due to the spin polarization transition discussed in the previous section 
for ID arrays. It can also be seen that the transition from the 6-electron to 7-electron ground state is forbidden 
at B = 1.3T and therefore the 7th peak is absent in Figs. |ll|(b) and |l2|(b). This is an example of the so-called 
"spin-blockade''_phenomenonJl3 The 8th peak at B = in the array is present in the conductance and partial width 
traces in Figs. |ll](a) and |l2|(a), although it is suppressed by a factor of approximately 100. Finally, for N = 1 the 
charge density on the corner dots is 10 times smaller at B — 1.3T, leading to the suppression of the first peak at 
B = 1.3T from its B = value in Fig. 0. 

We conclude that the spin polarization transition discussed in Sec. V for the case of the ID array leads to a similar 
suppression of the linear conductance peak heights in 2D arrays of quantum dots. 



VII. CONCLUSIONS 



We have shown that the formation of artificial molecular bonds due to interdot superexchange can drastically 
modify the low-temperature transport through coupled quantum dots. The resulting interdot spin-density- wave 
(SDW) correlations are strongly pinned by magnetic disorder, leading to a suppression of transport. These SDW 
correlations are destroyed in an applied magnetic field large enough to polarize all the electron spins, leading to a 
marked increase of the conductance at the spin-polarization transition (SPT). For a double quantum dot, this leads 
to a magnetoresistance proportional to (g/j,Bh/e 2 )U 3 /t 4 , where U is the charging energy of a quantum dot and t is 
the interdot hopping matrix element. Since U is typically at least an order of magnitude greater than t, we have 
termed this effect giant magnetoresistance. For larger ID arrays of quantum dots, the magnetoresistance was found to 
be proportional to (g ush / 'e 2 )U 2A r_1 / / t 2N at the JV-electron resonance when N is even, while saturating to a smaller, 
iV-independent value for N odd. This [/-dependence reflects the probability of an electron to tunnel all the way 
through the system while leaving the pinned, Neel ordered spin configuration of the ground state undisturbed. The 
giant magnetoresistance effect proposed here for coupled quantum dots is expected to be quite generic in quasi-one- 
dimensional systems with magnetic disorder. 

In addition to the giant magnetoresistance effect predicted for ID arrays of quantum dots with magnetic disorder, 
the SPT was shown to lead to large magnetoresistance effects in non-disordered ID and 2D systems and in systems 
with charge disorder. These many-body effects are likely to provide a fruitful area of research in the next generation 
of coupled semiconductor quantum dot systems. 
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FIG. 1. (a) Schematic diagram of a linear array of quantum dots, (b) A quantum dot array embedded in an Aharanov-Bohm 
ring, formed by connecting the left and right reservoirs in (a). 



12 



6 

—i 

E ' 

— 4 

in 

X 

2 



a n = 2 / 


/3 






□ 


□ n= 1 






A 












□ 


A 


A 


A 


□ 




□ 




□ 







0.0 0.5 1.0 1.5 

B (T) 

FIG. 2. Spin susceptibility \s = h~ 1 AS/AB at T — vs. magnetic field B for linear arrays of GaAs quantum dots with 
e 2 /C g = lmeV, C/C g = 0.5, and t — .05meV. Squares: 10 electrons on 10 dots (B c 1.5T); triangles: 8 electrons on 12 dots 
(B c » IT). 
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FIG. 3. The total wavefunction overlap r^/ r(i) = X^i E CT I (°n \4J0n-i) j 2 as a function of N for an array of ten quantum 
dots with pure Hubbard interactions U = Wt. The x coordinate of the iVth peak indicates the value of the polarization charge 
Q g — CgVg at which the iVth electron is added to the system. The suppression of with increasing iV is evident in the lower 
Hubbard band Q g /e < 1. The symmetry about Q g /e — 1 follows from the particle-hole symmetry of Eq. (|l|) 
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FIG. 4. Resonant conductance at T = in units of (e 2 //i)/(l + r' l '/r) as a function of the Zeeman splitting on dot 1 for 
7 _1 = 0, 3, 10, 30 (top to bottom). If an additional uniform external field B > J is applied, the conductance is restored to the 
value for 7 _1 = 0, leading to giant magnetoresistance. 
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FIG. 5. Persistent current through a double quantum dot embedded in an Aharanov-Bohm ring as a function of the gate 
voltage. Here t = 1, U = 10, U12 = 0, and T4i = Vfejv d = V. The current is expressed in units of J» = eV/2 1 ^ 2 h, the value 
on resonance in the absence of correlations and asymmetry. Solid curve: V = 1/32; dotted curve: V = 1/8; dashed curve: 
V = 1/2. 
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FIG. 6. Conductance vs. chemical potential /i through a linear array of 10 GaAs quantum dots with one spin-1/2 orbital 
per dot. e 2 /C g = lmeV, C = 0, t = O.lmeV, and T = 35mK. Splitting of the two Coulomb blockade peaks into minibands is 
driven by t. The suppression of the 5th peak in (b) is the result of a density-dependent SPT. 
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FIG. 7. Conductance vs. chemical potential p through a linear array of 6 GaAs quantum dots with one spin-1/2 orbital per 
dot. e 2 /C g — lmeV, C/C g = 0.5, i = .05meV, T = 120mK. Disorder 5t/i ~ 1 (tjf = tn) is present in the hopping matrix 
elements. The splitting of the Coulomb blockade peaks into multiplets is dominated by C; however, the effect of B is similar 
to that in Fig. 1. 
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FIG. 8. Conductance vs. chemical potential p through a linear array of 6 GaAs quantum dots with one spin-1/2 orbital per 
dot. e 2 /C g — lmeV, C/C g = 0.5, i = .05meV, T = 120mK. Spin-dependent disorder St/i ~ 1 (t,f =^ ti±) is included in the 
hopping matrix elements. Solid curve: 5 = 0; dotted curve: B — 1.3T. At 1.3T, the 2nd conductance peak is enhanced by a 
factor of 1600. 
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FIG. 9. Conductance vs. chemical potential through a linear array of 5 GaAs quantum dots with 5 spin-1/2 orbitals per dot. 
e 2 /C g = lmeV, A = 0.2 meV, T = .29K, C (n) /C g = 2"" 1 , and t„ = 0.05 meV(1.05) n (n = 0, . . . ,4). The energy gap between 
Hubbard minibands is not resolved for p > 9 meV (breakdown of CCB) . Note the quenching of magnetoconductance effects in 
the ballistic regime. 
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FIG. 10. Schematic diagram of a 3 x 3 array of quantum dots. 
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FIG. 11. The partial width = r^ CT |{Oiv I^L |Oat-i > | 2 as a function of N for the 3x3 array of quantum dots with pure 
Hubbard interactions U = lOt. The widths are plotted normalized by the partial width T% in (a) and F% in (b). The x 
coordinate of the iVth peak is given by E% — E%_ 1 . Each peak is labeled by N. The peak structure is discussed in the text. 
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FIG. 12. Conductance vs. chemical potential [i through the 3x3 array of quantum dots. The same parameters used as in 
Fig. 6. The peak structure is discussed in the text. 
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